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Abstract 

o 



We consider group-covariant positive operator valued measures (POVMs) on a finite 
dimensional quantum system. Following Ncumark's theorem a POVM can be imple- 
mented by an orthogonal measurement on a larger system. Accordingly, our goal is to 
find an implementation of a given group-covariant POVM by a quantum circuit using 
its symmetry. Based on representation theory of the symmetry group we develop a 
general approach for the implementation of group-covariant POVMs which consist of 
rank-one operators. The construction relies on a method to decompose matrices that 
intertwine two representations of a finite group. We give several examples for which 
the resulting quantum circuits are efficient. In particular, we obtain efficient quantum 
circuits for a class of POVMs generated by Weyl-Heisenberg groups. These circuits 
allow to implement an approximative simultaneous measurement of the position and 
crystal momentum of a particle moving on a cyclic chain. 
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1 Introduction 

General measurements of quantum systems are described by positive operator-valued mea- 
sures (POVMs) [1,2]. For several optimality criteria the use of POVMs can be advantageous 
as compared to projector valued measurements. This is true, e. g., for the mean square er- 
ror, the minimum probability of error [3], and the mutual information [4]. POVMs are more 
flexible than orthogonal von Neumann measurements and can consist of finite as well as of 
an infinite number of elements. An example for the latter is given in [5] where a POVM for 
measuring the spin direction is proposed. Here we restrict our attention to the finite case 
where a POVM is described by a set of positive operators which sum up to the identity. 
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Such a POVM is called group-covariant if the set is invariant under the action of a group. 
The example of POVMs for the Weyl-Heisenberg groups as well as the example in [5] show 
that POVMs are needed to describe phenomenologically the mesoscopic scale of quantum 
systems. They allow approximatively simultaneous measurements of quantum observables 
which are actually incompatible. For instance, the classical phase space of a particle can 
be approximatively reproduced by simultaneous measurements of momentum and position. 
Descriptions of quantum particles which have strong analogy to the classical phase space are 
helpful to understand the relations between the classical and the quantum world [6]. Also 
for several other tasks in quantum information processing the implementation of POVMs 
is of interest [7-9]. 

Neumark's theorem [10, 11] states that in principle every POVM can be implemented 
by an orthogonal measurement of the joint system consisting of the system and an ancilla 
system. However, the orthogonal measurement required by this construction may not be a 
"natural" observable of the joint system. One may need an additional unitary transform 
to obtain a reduction to a more natural observable which henceforth will be called the 
measurement in the computational basis of the quantum system. 

Therefore, the question arises how to actually implement a POVM in terms of a quantum 
circuit which itself is composed of a sequence of elementary quantum gates [12]. So far, 
only little is known about the implementation of POVMs even in quantum systems with 
a small number of dimensions. While some rather specific single-qubit measurements have 
been studied [4, 13, 14], not much is known about the general problem of how to implement 
a POVM by a unitary transform on the quantum register of a possibly larger space followed 
by an orthogonal measurement in the computational basis. 

When studying quantum circuits for families of POVMs questions about the complexity 
of the required unitary transforms arise. In some cases we can exploit the fact that they 
admit some additional symmetry. This leads to the study of group-covariant POVMs which 
has been studied extensively in the literature [4,15-17]. As a recent example we mention 
the construction of symmetric informationally complete POVMs by means of suitable finite 
symmetry groups [18]. 

The main contribution of this paper is a general method which computes an embedding 
of group-covariant POVMs into orthogonal measurements on a larger Hilbert space. A 
particular feature of the computed embedding is that it uses the symmetry. This in turn 
allows to apply known techniques for decomposing matrices with symmetry to the unitary 
matrices obtained by this embedding. For several cases this leads to families of efficient 
quantum circuits implementing the given POVMs. 

Outline. In Section |2] we briefly recall the definition of POVMs. In Section |3] we con- 
sider the decomposition of matrices that have a symmetry with respect to a group. This 
type of decomposition is a basic tool for our constructions. We also define group-covariance 
of POVMs with respect to a symmetry group and a group representation. Furthermore, 
we explain how POVMs with this group-covariance are related to so-called monomial rep- 
resentations of the symmetry group. In Section 0] we explain the general scheme for the 
construction of a unitary transform that implements a group-covariant POVM. The basis for 
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this construction is the analysis of the intertwining space between the group representation 
that is given by the group-covariance of the POVM and the monomial representation. This 
is the starting point for methods using fast quantum Fourier transforms as described in Sec- 
tional Finally, in SectionHOwe give several examples of implementations of group-covariant 



Notations. We denote the field of complex numbers by C. The group of invertible rax ra 
matrices is denoted by GL n (C) and the subgroup consisting the unitary n x n matrices is 
denoted by U{n). We denote the identity matrix in U(n) by l n . If not denoted otherwise 
all matrices are matrices over the complex numbers. The cyclic group of order n is denoted 
by Z n . Representations are denoted by small Greek letters, e.g., tp, ip etc. By abuse of 
notation we also denote the trivial representation of degree n (i.e. dimension ra) by l n . The 
base change of a matrix A with respect to a matrix B is denoted by A B = BAB^ . The 
direct sum of matrices and representations is denoted by A © B and <p®ijj and the tensor 
product is denoted by A © B and <p (g> if), respectively. We make frequent use of the Pauli 
matrices 



A diagonal matrix with diagonal entries Ai, . . . , A n is abbreviated by diag(Ai, . . . , A n ). We 
denote the symmetric group on ra symbols by S n . To each permutation a G S n naturally 
corresponds the permutation matrix Y2i l°"00)(*l- By abuse of notation we identify a with 
the corresponding permutation matrix. We often use the permutation matrix S m which 
corresponds to the m-cycle (1,2, ... ,m) and the matrix T m = diag(l, u m , . . . which 
contains the eigenvalues of S m . The basis states of an n-qubit system correspond to binary 
strings of length n. Quantum circuits are written from the left to the right, and the qubits 
are arranged such that the most significant qubit (characterizing the left-most symbol of a 
binary string) is on top. Throughout the paper a matrix entry "•" stands for zero. 

2 POVMs and orthogonal measurements 

A POVM for a quantum system with Hilbert space C d is a set P = {A\, . . . , A n } C C dxd of 
non-negative operators, where Ylk = Id- For a more general definition for POVMs with 
an infinite number of operators we refer to [19] . For example, the set of matrices 



where uj = exp(27ri/3) is a third root of unity, defines a POVM on a system with corre- 
sponding Hilbert space C 2 . Suppose that the state of the system is described by the density 
matrix p G C dxd . Then for a general POVM the probability pk for the result k is given by 
Pk = tr(pAk). An orthogonal measurement is a POVM with mutually orthogonal operators 
Ak, i- e., we have that A^Ai = A[Ak = for k ^ I. 



POVMs. 
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In the following we restrict ourselves to rank-one operators A k = \^ k ){^ k \. Note 
that the POVM vectors \^k) need not be normalized and that the restriction to opera- 
tors of rank one is for some applications justified by Davies' theorem [15]. It states that 
we can always find a POVM with rank-one operators that maximizes the mutual infor- 
mation. The example P2, which consists of three rank-one operators, can be written as 
P 2 = |* 2 >(*2|, |*3><*3|}, where 

are the corresponding POVM vectors in C 2 . Neumark's theorem [11] states that it is 
possible to implement a POVM by reducing it to an orthogonal measurement on a larger 
system. We briefly recall this construction. Let P = {A k } = {\ty k ) k \} be a POVM 
with n operators that acts on the Hilbert space C d . For n > d the vectors \^ k ) cannot 
be mutually orthogonal. Consequently, we have to extend the system by at least n — d 
dimensions in order to define an orthogonal measurement with n different measurement 
outcomes. We want to implement an orthogonal measurement P = {A k } = {\^k)(^k\} on 
the system with n dimensions such that P corresponds to the POVM P on the subsystem 
with d dimensions, i.e., p k = tr(pA k ) = tr(pA k ). Here p = p © n -d £ C nxn where n -d 
denotes the zero matrix of size n — d is the embedding of the state into the larger system. 

We write the POVM vectors as columns of the matrix M = . . . G C dxn . 

In the following we refer to M as the defining matrix for the POVM P. Now, the operators 
A k = \y k ){$ k \ G C nxn with = © l^fc) are the columns of the matrix 

Note that M can be an arbitrary unitary matrix which contains M as upper part of size 
dx n. Since Pisa POVM we have MM 1 " = J2k \^k)(^k\ = Y.k A k = Id, i-e., finding a 
suitable M is always possible. For example in case of P2 we obtain the defining matrix 



M = -1= ( \ \ 1 ) G C 2x3 

and one possible choice for M is to add the row given by (l/v / 3)(l,^,^ 2 ). Hence the 
rank-one projectors corresponding to the orthogonal measurement M are 




*x) = -= 1 , |* 2 ) = -= uj 2 ,and |*3) 




The probability distribution p k = tr (pA k ) of the constructed orthogonal measurement 
equals the distribution p k of the original POVM since 

ft = tr (pA t ) = tr ((, ffi ( W<J;J W<Jjj )) = tr (M„) = P, 
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The embedding into a larger system can be realized by using an ancilla register of a quantum 
computer. It consists of I qubits such that 2 l > n — d. They are initially in the state |0 . . . 0). 
Then the space C d <g> |0 . . . 0) is the subspace where the POVM acts on and C d (C 2 )® 1 
is the extension. The density operator p acts on an n dimensional subspace of the joint 
system consisting of the original system and the ancilla register. In the following we will 
assume that also the system space C d is embedded into the state space of some qubits. 

As explained above, we can implement the POVM with corresponding matrix M by 
applying the unitary transform AF to the initial state p of the joint system followed by a 
measurement in the computational basis. Note that for the special case where the columns of 
M are already orthogonal we have that M = M. In this case by implementing the matrix 
Mt followed by a measurement in the computational basis we can perfectly distinguish 
between the columns of M. 

In principle, the construction of an appropriate matrix M is simple since we just have 
to find mutually orthogonal rows that lead to a unitary matrix. However, k qubits allow 
POVMs with n = 2 k operators. Hence the size of M is exponential in k. The complexity to 
implement a unitary matrix on k qubits can be upper bounded by 0(4 fc ) [20] and a generic 
element of lA(2 k ) will indeed require an exponential number of elementary transforms (e.g. 
one- and two-qubit-gates). Therefore we are interested in the construction of a matrix M 
that can be implemented efficiently, if such a construction exists at all. While finding effi- 
cient factorizations is a hard problem in general, the situation becomes easier in some cases 
where we are given the additional structure of a group-covariant POVM. In the following 
sections we will give a definition of group-covariance and the related notion of symmetry. 
Later, we exploit the symmetry of the matrix M and give several examples of POVMs that 
have efficient quantum circuit implementations. 



3 Group-covariant POVMs and matrices with symmetry 

In the following we give a precise mathematical definition of the notion of symmetry of a 
matrix M G C mxn . Later we define group-covariance of a POVM and show that the group- 
covariance in a natural way leads to matrices with symmetry. For the necessary background 
on finite groups and representations we refer to standard textbooks such as [21,22]. 

We start with a finite group G and a pair (ip, ip) of matrix representations of G which 
are compatible with the size of M, i. e., ip : G — > GL m (C) and tp : G — > GL n (C). Following, 
[23,24] we call the triple (G,ip,ip) a symmetry of M if the identity f(g)M = Mip(g) holds 
for all g G G. Sometimes we abbreviate this by using the shorthand notation ipM = Mip. 
Note that if M is not a square matrix the representations <p and ip have different degrees. 

To give an example we let oj = exp(27ri/3) and let a, (3, 7 G C. Then for all j G {0, 1, 2} 
we have that 

(1 • \ J / a a a \ (a a a \ / • ■ 1 

■ u ■ i L- U :I 111 • 

u) 2 J \ 7 ju> 2 juj J \ 7 7CJ 2 juj J \ ■ 1 • 
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Hence we obtain a symmetry which is given by the cyclic group Z3 = {0, 1, 2} together with 
the two representations tp, a : Z3 — > 14(3) given by tp(l) = diag(l, u>, u> 2 ) and a(l) = (1, 3, 2). 

Note that given two representations <p, tp of a group G the set of all matrices M which 
fulfill tp(g)M = Mip(g) for all g G G is a vector space. It turns out that the matrices in 
this vector space have a special form. Hence we explore its structure in more detail in the 
following. 

Definition 1 (Intertwining space) Let G be a group and let tp, ip be representations of 
G of degrees n and m, respectively. Then 

Int(v9,V0 := {M : tp(g)M = Mip(g), for all g G G] 

with M G C nxm is called the intertwining space of tp and ip. 

In the following we denote by tpi,...,tpk a list of all pairwise inequivalent irreducible 
representations of G. Recall that for any representation of a finite group it is always 
possible to find a base change such that the corresponding representation is a direct sum of 
irreducible representations [22]. For representations which are completely decomposed into 
a direct sum of irreducibles the structure of the intertwining space is known. This is the 
content of the following theorem which follows directly from Schur's Lemma (see [25, Section 
§29]). 

Theorem 2 Let G be a finite group and tp = 0f =1 (ln i ® </>i) and ip = 0f =1 (l 

two representations of G which have been completely decomposed into pairwise inequivalent 

representations (fi, i = 1, . . . , k. Then the intertwining space of if and tp has the following 

structure: 

1nt(<p,1>) = (C" lXmi ® i deg(¥Jl) ) © . . . e (C n ** m * ® i degiVk) ). 

A matrix A is called block permuted if there are permutation matrices P and Q such 
that PAQ = B\ © . . . © Bk, where B\, . . . , Bf. are (rectangular) matrices. For all n, m, k G N 
there exist permutation matrices P nt m,k an d Q n ,m,k such that for all A G C nxm we have 
Pnmk(A® lfc)Qnmfe = lfc ©^4- Hence we have shown that the elements of the intertwining 
space of completely reduced representations are block permuted. 

We continue with an easy observation which turns out to be essential for the approach 
of extending the symmetry of a given group-covariant POVM to a measurement on a larger 
space. Suppose that M G Int(y, tp) and that the matrices U and W decompose the represen- 
tations tp and tp into the direct sums, i. e., UtpU^ = <pi® . . .®<p n and VipV^ = ipi, © . . .®ipm- 
Then we can rewrite tpM = Ma as 

© . . . © tp n )UM = MW\ipi © ... © tp m )W. 

Multiplying this from the left by U and from the right by W' shows that C := UMW^ is 
an element of the intertwining space Int((/?i © . . . © tp n , ipi © . . . © ip m ) of two completely 
reduced representations. In particular, we can apply Theorem|5]to determine the structure 
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of C . In particular we obtain that C is block permuted and the size of the blocks depend 
on the multiplicities and degrees of the irreducible representations contained in ip and ip. 

Matrices with symmetry arise naturally in context of group-covariant POVMs. We first 
give a definition of these POVMs and then establish a connection between the notions of 
group-covariance and symmetry. 

Definition 3 (Group-covariant POVMs) A POVM P = {A 1 ,...,A n } C C dxd with 
A k 7^ A\ for k 7^ / is group-covariant with respect to the group G if there exists a projective 
unitary representation ip : G ^U(d) with ip(g) A k ip(g)^ G P for all g G G and all k. 

Note that a group-covariant POVM is also group-covariant for all subgroups H < G and 
the restriction of the representation ip to H. As a special case, the choice of the trivial 
subgroup H = {1} means that we do not use the symmetry of the POVM at all. 

A minor complication arises due to the fact that while the notion of symmetry of matri- 
ces relies on ordinary, i.e., non-projective representations, the definition of group-covariant 
POVMs relies on projective representations. Therefore, we need a construction which al- 
lows to transform the projective representation of the symmetry group of a group-covariant 
POVM into a non-projective representation. This connection is established using so-called 
central extensions which is a method going back to I. Schur. We briefly recall this construc- 
tion (see also [22, Lemma (11.16)]). Let p : G — > GLrf(C) be a projective representation 
of the group G. More precisely, we have <p(gh) = r y g h ( f(9) < p(h) for g,h G G, where 7^ is a 
factor system. Let H = (7^ : g, h G G) be the group generated by the j g h- We consider 
the group G consisting of the elements (g, h) with g G G and h G H. The multiplication 
of two elements (g, h) and (</, h!) of G is defined by (g, h)(g', h') = (gg' ,j gg ihh'). Then the 
map (p((g, h)) = hp>(g) is a representation with <p~((g, 1)) = ip(g), i. e., the representation <p 
equals ip on the elements (g, 1) and the group G is a central extension of the group G. 

In the following we always assume <p to be a non-projective representation of the sym- 
metry group G by this construction. This is justified since the set of POVM operators does 
not change by switching from G to a central extension G because scalar multiples of the 
identity operate trivial under conjugation. 

We now analyze the structure of the matrix M corresponding to the group-covariant 
POVM P = {\^k)(^k\} with rank-one operators. Note that the phases of the vectors 
can be chosen arbitrarily without changing the POVM. Let ip : G —> U(d) be the 
representation corresponding to the symmetry of P. We then have the equation 

p{g)\^ k )^ k \ip( g y = \y <g)k ){y <g)k \ 

where n : G — > S n denotes a permutation representation of the group G. Indeed, the 
equation l^^-X*^! = \^7T( g )k)(^n( g )k\ implies = |*fc)(*fc| by conjugation 

with tp(g)^ since Aj 7^ A k for j 7^ k. Therefore, the map ir(g) is injective for all g G G. 
Since an injective map on a finite set is also surjective the map ir(g) defines a permutation. 

Next, we consider the action of <p on the columns of the matrix M. As stated above the 
columns \^ k ) of M can have arbitrary phase factors. The action of ip(g) on the columns 
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of M can be described by the equation <p(g)\tyk) = ei ^ 9 ' k ^\^ir(g)k) where 4>(g,k) depends 
on k, g and the fixed phase factors of the vectors \^k)- We identify the columns \^k) with 
a basis bk of the vector space C n in order to construct a representation that describes the 
action of </? on the columns of M. With this identification the action of (p{g) corresponds 
to the map bk >—>■ e^^'^b^^k- 

By writing down the matrix corresponding to this map, we see that in each row and 
each column there is precisely one entry different from zero. Matrices having a structure 
like this are called monomial matrices 1 [25, Section §43]. Whenever the images under a 
representation consist entirely of monomial matrices, we denote this with an underscript, 
i.e., we write <^mon(<7)- Now, the two representations <p and <p m0 n define the symmetry 
ipM = Mip 

mon °f the matrix M . The monomial representation <p mori acts on the columns 
of M. For each g G G it permutes the columns of M and multiplies each column with a 
phase factor. 



Example 4 As an example in two dimensions we consider the following POVM: 



P = 



\a\ 

a(3 



af3 



a 



-a{3 



-a(3 

«|2 



|/3|_ 2 <*0_ 
a(3 



a 



—a(3 



-a(3 



cc 



2x2 



with a, P G C and \a\ 2 + \(3\ 2 = 1/2. Then P is covariant with respect to Z2 x Z2. The 
corresponding projective representation ip : 7Li x — ► U{2) is defined by the equations 

99(0, 0) = 1 2 , y(0, 1) = cr 2 , 99(1, 0) = a x , tp(l, 1) = a z a x 

where (0, 0), (0, 1), (1, 0) and (1, 1) denote the elements of the group Z 2 x Z 2 . 

For this projective representation of Z 2 x Z 2 a simple computation shows that the 
central extension G of Z 2 x Z 2 is isomorphic to the dihedral group with eight elements. In 
the following it is sufficient to consider the definition of the representation on the elements 
((0, 1), 1) and ((1, 0), 1) since these elements generate G = {(g, h) : g G Z 2 x Z 2 , h G {±1}}. 
We can choose 

a a (3 f3 \ ^ c2x4 
(3 —[3 a —a 



M 



or a matrix with the same columns (up to an arbitrary phase factor for each column). This 
leads to a symmetry group given by the monomial representation 



¥>mon((0, 1), 1) 



/ • 1 
1 • 

V • • 



\ 



/ 



and <^mon((l,0), 1) 



V 



1 



-1 



, . The modified 



For a different choice of phase factors we obtain another representation ip B 
pair of representations tp, ip mon also defines a symmetry of M. 

1 Note that this terminology is somewhat unfortunate since it has nothing to do with the monomials of 
which a polynomial is comprised of. Still it is the standard terminology used in representation theory. 
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An important special case of group-covariant POVMs are group- generated POVMs which 
we describe next. Let G be a group and p : G — > <C dxd an (ordinary) unitary representation. 
A group-generated POVM is described by the POVM vectors <p(g)\ty) for g E G and an 
initial vector \^>) E C d . The corresponding operators of the POVM are given by A g = 
p{g)\^>) {^\p(g)^ for g G G. In other words, all POVM vectors are obtained by the initial 
vector under the operation of the group G, i.e., they form an orbit. Obviously, a 
group-generated POVM is a group-covariant POVM with a single orbit under the action 
of the group. With this construction, the phase factors of the POVM vectors ip(g)\^f) 
are fixed by the phase factor of the initial vector |*). The phase factors e 1 ^ 9 '^ of the 
monomial representation (p mon corresponding to <p equal 1. As a consequence, the monomial 
representation p mon equals the regular representation of G where we have to consider a fixed 
order of the elements of G. 

Note that the operators {(p(g)\^) (^\p(g)^ } in general do not define a POVM for arbi- 
trary representations <p and initial vectors |\&). However, if <p acts irreducibly one has (after 
appropriate normalization) for every vector the equation ^2 ge Q p(g)\^} (^\p{g)^ = 1<2- 



4 Construction of the orthogonal measurement 

Following the previous section we can arrange the vectors which correspond to the elements 
of a POVM with rank one projectors into the columns of a matrix M. We have seen 
that in case of a group-covariant POVM the matrix M G C dxn always has the symmetry 
(pM = Mip mon where ip is the given representation and ip mon is a monomial representation. 
Both representations are representations of the symmetry group of the group-covariant 
POVM. We know that both representations are equivalent to direct sums of irreducible 
representations. Hence we can find unitary matrices U and W such that UpU^ = ipi © . . . © 
ip n and Wp mon W^ = <7i©. . .©cr m where the pk and the 07 denote irreducible representations 
of the group G. In general, we can write the equation (pM = Mip mon as 

U\<p! © . . . © (f n )UM = MW ] {ai © ... © a m )W. 

This is equivalent to C = UMW^ £ T := Int(^i ©...©</?„, 01 © ... © a m ). Conversely, 
a matrix C which is contained in this intertwining space and has orthogonal rows defines 
(up to an appropriate normalization) a group-covariant POVM with corresponding matrix 
M = U^CW. 

For a given matrix M £ <C dxn we now consider the construction of a unitary matrix 
M G U(n) such that M contains M as upper part, i. e., we are looking for a matrix M such 
that 

' M 



M 



N 



where TV G c( n ~ d ) xn , \ n addition to this we intend to get the symmetry (</?©<//)M = Mp mon 
with an appropriate representation p' : G — ► U(n — d). If we succeed in constructing an 
appropriate representation ip' and matrix M then we have the equation p(£>p' = Mp mon M\ 
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i. e., the representation ip®ip' has to be equivalent to y) moil . In other words, each irreducible 
representation of G is contained the same number of times in <p ffi tp' and in ip mon . Further- 
more, from the decompositions (U © \ n -d)(<p © tp')(U^ © l n -d) = °V(i) © • • • © <?V(m) and 
W(p m0 nW* = o"i © ... © o~ m we obtain that 

(U © 1)MW ] £ f := Int(V T( i) © ... © a r(m) ,ai © . . . © a m ) C C" x ". (1) 

The permutation r used in eq. Q is a suitable reordering of the irreducible representations. 
The structure of the intertwining space T is known from Theorem [2] since we can compute 
the irreducible representations aj from (p m0 n- 

In the following discussion we consider the construction of (p' and M. Our goal is to show 
that the construction of ip f that makes U(pU^ © <p' equal to Wip mon up to a permutation 
r of the irreducible components is always possible. 

Important for the extension of M to M will be the following theorem which characterizes 
the relations of two representations in case there is an intertwiner of maximal possible rank. 
Recall that ipi is a constituent of ip2 if and only if there is a base change U such that 
U Tp2(g)U = tpi(g) © tp'i(g) where ?/4 is a representation of G. 

Theorem 5 Let G be a finite group and let ipi,tp2 be representations of G of degrees d\ = 
deg(V ; i) and c?2 = deg(?^2), respectively. Let M £ C lX 2 be a matrix with ifji(g)M = 
Mip2(g) for all g £ G and rk(M) = deg(-0i). Then ip\ is a constituent of tp2- 

Proof: Let M be such that il)\{g)M = Mip2(g) and let (pi,...,tpk be a complete set of 
pairwise inequivalent irreducible representations of G. Since ^i, ip2 are representations of 
a finite group over the field of complex numbers we find unitary matrices U, W such that 
UipiU^ = ©j =1 mnpi and WfaW^ = ©f =1 ?ij9?j, where the multiplicities m« and n« are 
non- negative integers. We have to show that actually mj < n» for all z = 1, . . . , k. 

From ipiM = Mip2 and by the choice of U and W we obtain that m,iipi)(U MW>) = 
(UMW^)(®nnpi), i.e., we have that UMW* £ lnt(0 t fc =1 mm, 0*^ n Wi)- % the re- 
marks following Theorem|2]we know that there are permutation matrices P and Q such that 
M := P(UMW^)Q = (l deg{vi) ®B 1 )(B...(B(l degivk) ®B k ) where each B< £ C m > xn *. Mul- 
tiplication with invertible matrices preserves the property that M and hence also Mq have 
full rank (given by deg(V>i))- On the other hand we know that the rank of a block diagonal 
matrix is given by the sum of the ranks of the blocks. Hence rk(Mo) = Y2i=i deg(ipi) -rk(Sj) 
which shows that each B{ must have full rank. Since Bi is an mjXnj matrix this in particular 
implies that mi < n^. This shows that ipi is a constituent of ip2- n 

We now use Equation Q to construct the matrix M for the implementation of a 
group-covariant POVM. Having determined U and W we can compute the matrix C = 
UMW^ £ lnt(UipU\ WifmonW^). The number of times each irreducible representation 
has to occur in <p' can be computed. Since the structure of the intertwining space T = 
In{,{UipU^ © tp', W(p mon W^) is known we can extend C to an arbitrary unitary matrix C 
of the intertwining space T. This extension is always possible since both representations 
U<pW © tp' and Wip 

mcmW contain each irreducible representation the same number of 
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times. The matrix C defines some of the rows of A. Since M defines a POVM the rows 
are mutually orthogonal. Consequently, the matrix components of C corresponding to an 
irreducible representation can be chosen under the constraint that they are orthogonal. We 
now have that for any V £U(n — d) the matrix M = (U^ V^)CW yields a unitary that 
extends the matrix M and has the symmetry we wanted to construct. 

Hence, we obtain the following algorithm to construct an orthogonal measurement which 
realizes the given POVM and preserves the symmetry. 

Algorithm 6 Let P = {A-y, A n } C C dxd be a POVM. Then the following steps imple- 
ment P by a von Neumann measurement on a larger space. 

1. Write the rank-one operators A^ = ^^(^j of the POVM as columns of the matrix 

M G C dxn . 

2. Determine an appropriate symmetry group with corresponding representation ip : 
G^U{d). 

3. Compute the monomial representation <p mon : G — ► U(n). 

4. Find a matrix U G U{d) that decomposes ip into irreducible representations where 
equivalent ones are equal. 

5. Find a matrix W 6l/(n) that decomposes f mon into irreducible representations where 
equivalent ones are equal. 

6. Construct the representation <p' such that UpU^ © ip' is equal to Wip mon W^ up to a 
permutation r of the irreducibles. 

7. Construct C G U{n) that contains C = UMW^ G C dxn as upper part and is in the 
intertwining space T of UtpU^ © <p' and VF^mon^^- 

8. Choose an arbitrary unitary matrix V £U(n — d). 

9. Compute M = (W © W)(7VF G U{n). 

Then implements the POVM P by a von Neumann measurement on a larger space, i. e., 
for any state p on the original ci-dimensional system we have that pk = tr(pAk) = (^>k\p\^k)- 
Here denote the rows of M and p = p © n _d is the embedding of p to a state of an 
n-dimensional system. 

Example 7 We consider the example of the previous section with the matrix 

M=(° " P MgC 2 * 4 
y p —p a —a J 

and the group G = {(g,h) : g G Z2 x Z2,/i G {±1}} which is isomorphic to the dihedral 
group of order eight. The representation (p : G — > U(2) is given by <^((0, 1),1) = a z 
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and <£>((1,0), 1) = a x . We have U = I2 and UpU^ = tp since the representation </? is 
already irreducible. An elementary computation shows that the corresponding monomial 
representation p nion is given by 



WV mon ((0,l),l)^t = 
with the unitary matrix 



/ 1 • 

• -1 

• 1 

v • • 

V2 



■ \ 

"I / 
/ 1 

V 1 



and Wp moQ ({l,0),l)W^ 



\ 



1 



1 • -\ 

• 1 -1 

• 1 1 

1 • •/ 



€W(4). 



T = lnt((f (£>p,p(Bp) 



■ Xij G C > c 



-•4x4 



Therefore, p mon contains the irreducible representation p twice, i. e., W(p mon W^ = p © p>. 

With the matrices M G C 2x4 ,[7 G W(2), and G W(4) as above we find that C = 
UMW^ = y/2 (a (3 ) (8> I2 G C 2x4 , which is an element of the intertwining space 

Int(<£, W^monW^) = Int(y>, p © 93) . 

Since we have W^ mon W^ = 99 © 99, we have to choose p>' = p. The intertwining space T is 
given by 

/ An • A12 • \ 
An • A12 
A21 • A22 

\ • A21 • A22 / 

In our example, the matrix C = UMW^ defines the first two rows of the matrix C G T = 
lnt(p © p,p © p). 

In particular, we have the equations An = y/2a and A12 = V%P- For example, it is 
possible to choose A21 = \/2 (3 and A22 = — V^u for a, (3 G C to obtain the unitary matrix 

/ a ■ (3 \ 
■a (3 
(3 • —a 
\ ■ P ■ -a ) 

which has the symmetry (p © p)C = C{p © p). With M = (U^ © V^)CW and V = 1 2 we 
compute the matrix 



C= V2 



G W(4) 



M 



/ 


a 


a 


P 


P \ 




/? 




a 


—a 




/3 


P 


—a 


—a 


V 


—a 


a 


P 


~P J 



G W(4) 



that contains M as upper part and has the symmetry (93 © p)M = Mp B 
unitary matrices V G U{2) give rise to possible extensions M. 



Note that all 
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5 Efficient implementations of group-covariant POVMs 



From the general construction of a von Neumann measurement which realizes a given POVM 
using the symmetry of the POVM we now turn to the question of decomposing the unitary 
M into gates. This can be seen as a first step towards the more general question of how 
POVMs can be implemented efficiently on a quantum computer. 

When speaking about the efficiency, we mean the cost of implementing the POVM as a 
von Neumann measurement on a larger Hilbert space, i. e., the number of elementary gates 
we need to actually implement the necessary unitary operation on this bigger space. First 
note that the discussed construction of M has several degrees of freedom: 

• The matrix C that contains C as upper part can be chosen arbitrarily. The matrix C 
has to be a unitary matrix in the intertwining space T. 

• The matrix V &U(n — d) can be an arbitrary unitary matrix. 

• The order and phase factors of the POVM vectors in the matrix M can be chosen 
arbitrarily. However, it must be possible to deduce the applied POVM operator from 
the result of the orthogonal measurement efficiently. 

• The permutation r of the irreducible representations in U (fU^ © <p' can be chosen 
arbitrarily. 

• The symmetry group G can be restricted to subgroups H < G which might lead to 
different realizations of the POVM. 

The constructions depend on the symmetry group G we consider for the POVM. Some- 
times, we can obtain simple implementations by restricting the symmetry group to a sub- 
group H < G. If we consider a subgroup H of G and construct the POVM with respect 
to H we have several changes in the construction compared to the construction with the 
group G. On the one hand, the number of occurrences of the irreducible representations in 
Vmon increase. On the other hand the number of inequivalent irreducible representations 
of the symmetry group decreases. Consequently, the matrices of the intertwining spaces 
are more complex since there are more irreducible representations in ip and y? mon that arc 
equivalent. As a tradeoff we have that the complexity of the transform W decreases. The 
circuits constructed in [14] show that the restriction of the symmetry group to a cyclic 
subgroup can lead to efficient algorithms in some cases. 

Let G be a finite group and {<pi, . . . , ipk} a system of representatives for the irreducible 
representations of G. Let the coefficients of these representation be indexed by the list 
V := [(m;i,j),l < m < k, 1 < i,j < deg(<p m )]. Furthermore, let the elements of G be 
indexed by the list L. Then the matrix l/y/\G\ (\J deg(ip m ) <Pm(g)ij)(m;i,j),g is unitary and 
is called a Fourier transform (or DFT for short) for G [26,27] (with respect to L and L'). 

For several groups it is known how to realize a DFT efficiently on a quantum com- 
puter [28-30]. In these cases the symmetry (p mon can be decomposed efficiently whenever 
we have that (i) <£> mon is a regular representation of G and that (ii) the DFT for G can 
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be computed efficiently. Note that the computational complexity of this von Neumann 
measurement depends essentially on the complexity of implementing DFT^r in terms of 
elementary quantum gates. Hence we obtain several families of POVMs for which the 
monomial representation </? mon can be decomposed efficiently. The complexity of the corre- 
sponding POVM then depends on the remaining matrices C, U, and W used in Algorithm 

m 



6 Examples 

In this section we apply the methods discussed in the preceding sections to some examples 
of group-covariant POVMs. We exploit the symmetry of group-covariant POVMs with re- 
spect to cyclic groups, dihedral groups, and Weyl-Heisenberg groups in order to construct 
quantum circuits for the implementation of these POVMs. Quantum circuits for the im- 
plementation of group-covariant POVMs on a single qubit with respect to the cyclic and 
dihedral groups are also discussed in [14]. 



6.1 Cyclic groups 



Let Z n = {0, 1, . . . , n — 1} be a cyclic group with n elements and let u = exp(27ri/n) 
be a primitive nth root of unity. On a d-dimensional Hilbert space we consider a group- 
generated POVM with respect to the representation ip : Z n — ► ti{d) that is defined on the 
generator by </?(l) = diag(l, uo, u/ 
the elements (p(g)\^>) for g G 



vector |^) 
matrix 



lAMi, 



M 



in 



iy 

( i 
i 



V i 



With an appropriate initial vector |^) G C 
7 j n define a POVM. In the following, we only consider the 
G C d . This vector leads to the POVM with the defining 

\ 

(2) 



UJ 



1 



UJ 



uJ 



1 

n-1 



l-l 



2{d-l) 



UJ 



(n-l)(d-l) 



•*dxn 



The matrix M G C rfxn has the symmetry ipM = M(p mon where </? m on(l) = (1 5 2, . . . , n). 
The representation tp mon is the regular representation of the cyclic group where the elements 
are ordered as [0, 1, . . . , (n — 1)]. With the Fourier matrix 



1 



in 



UJ 1 



n— 1 



j,k=0 



eU(n) 



we can write F n (/? mon (l) Fn = diag(l, u, to 2 , . . . , w n_1 ). This shows that the Fourier trans- 
form decomposes the regular representation of 7L n into a direct sum of irreducible represen- 
tations. 

According to the preceding discussion (and notation) we have that U = Id and W = F n . 
As a consequence we have the equation C = UMW^ = MFn- More precisely, we have 
C = MF n = diag(l, 1, . . . , 1) G C dxn . 
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We now consider the construction of the matrices C and M. The representation <p : 
Z n —* lA{d) with 93(1) = diag(l, u, uj 2 , . . . , uj l ) contains the irreducible representations 
1 1— > {uj k ) for all k G {0, . . . d — 1}. The representation F n (p mon Fn : Z n — > £Y(n) with 
^nV5mon(l)^n = diag(l, w, u; 2 , . . . , w n_1 ) contains the irreducible representations 1 1— > 
for all A; G {0,1,..., n — 1}. Following Algorithm El from Section^ we choose cp' with 
y/(l) = diag(w rf , . . . ,uj n ~ l ) in order to obtain </?©(// = -F n (/? mon F^. Since each irreducible 
representation 1 1— > with k G {0, 1, ... ,n — 1} has dimension one and the irreducible 
representations defined by 1 1— ► (oj k ) are inequivalent for different k we have the intertwining 
space 

f = M(p © ^, F n93mon Ft) = {diag(Ai, . . . , A„) : \j G C} C C nxn . 

We have to find a matrix (5 G U(n) in the intertwining space T that has the matrix C G C dxn 
as upper part. As stated above, the matrix M G C dxn defines Xj = 1 for j G {0, 1, . . . , d— 1}. 
Since C has to be a unitary matrix we have to choose Xj with the absolute value \Xj\ = 1 
for j G {d, . . . , n — 1}. 

In order to simplify the matrices we set Xj = 1 for all j G {d, . . . ,n — 1}. With these 
elements Xj we have the equation C = l n . Furthermore, we choose V = l n -d hi Algorithm 
El from Section ^ leading to U © V = l n . Consequently, we obtain the equation 

Mt = W 1Q\U ®V) = Fll n l n = Fl 

This equation shows that the inverse Fourier transform = Fn is a unitary transform that 
implements the group-covariant POVM with defining matrix (J2J. Recall that for n = 2 k 
where k G N the Fourier transform can be implemented efficiently on a qubit register [31,32]. 

6.2 Dihedral groups 

Let = ( r ,s : r m = l,s 2 = l,srs _1 = r _1 ) be the dihedral group [33] with n = 

2m = 2 k+1 elements for a fixed m = 2 k > 4. The element r denotes the rotation and s the 
reflection of the dihedral group. We consider the irreducible representation ip : Dim ^(2) 
that is defined by 



<p(r) = ( u t ) and v (s) = ( J J ) 



The element uj = exp(2-7ri/m) is an mth root of unity. For a,f5 G C with |a| 2 + \(3\ 2 = 1/m 
we consider the POVM with the corresponding matrix 

M=( ° * P t-i)^ Xn - 

\ p ... poj a ... aoj m 1 ) 

The matrix M G C 2xn has the symmetry ipM = M<p mon where (p mon is defined by the 
equations (^ mon (r) = 1 2 © ooS m 2 and ^ m0 n('S) = <J X ® F^T m . The matrices S m ,T m G C mxm 
are defined by the equations (indices are taken modulo m) 

m—l m—1 



i=0 i=0 
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and F m denotes the discrete Fourier transform defined in the previous section. In order to 
decompose </? m0 n into irreducibles the following permutation Qk is useful. Denoting by x 
the binary complement of the binary vector x of length k we define Qk '■ \x, 0) i— > \x, 0) and 
Qk '■ \x, 1) i— > \x, 1). Furthermore, we introduce the representations <pi defined by 

W(r)=(o w °-«) and w (s) = ( J J 

With this notation we have ip = ipi. The two-dimensional representations y?j are irreducible 
and inequivalent [33] for different I G {1, . . . , m/2}. Now, using the base change W := 
Q m (l 2 <8> i4) G C nxn we obtain that 

where ip is a direct sum of all representations yjj with odd j. The first component of ip is 
, the other components ipj appear in a specific order which is irrelevant in the sequel. We 
choose the representation 

tp' = ip'®ip@ip®ip, 

where if/ is obtained from ip by dropping <p\. This leads to ip © ip' = W(p mon W^. The 
matrix C = MW ] = (^fma . . . 0|V^/?0 . . . 0) <g> 1 2 € C 2xn defines the first two rows of 
the intertwining matrix C we want to construct according to Algorithm |f3 from Section |1J 
A possible extension of the intertwining matrix C G C 2xn to a unitary matrix C G U(n) is 
C = A <8i l m /2 with the matrix 



-4 = ^(9 



According to Algorithm H3 from Section ^ we have to define the matrices U G IA{2) and 
1/ G U(n — 2). The equations <p = <p>\ and W(p mon W^ = ((p\ @ if)') @ if> @ if) @ if> show that 
U = 12- Furthermore, we choose V = l n -2- Then we have the matrix U © V = l n - To 
summarize, we have to implement the matrix 

Aft = W tc-t = ( l2 g, F m )Q fc (,4 t U) G W(n) 

in order to measure the POVM corresponding to the dihedral group D m . The scheme of 
the circuit corresponding to Aft is shown in Figure ^ 



6.3 Weyl-Heisenberg groups 

In the following we introduce the finite Weyl-Heisenberg groups which are matrix groups 
acting on a finite dimensional vector space. For our purposes we consider vector spaces of 
dimension m = 2 k only, where k > 2. Then the Weyl-Heisenberg group G m is the group 
generated by the matrices S m = (1,2, ... ,m) and T m = diag(l, u, uj 2 , . . . , uj m ~ l ) where 
ijj = exp(27ri/m) G C is a primitive mth root of unity. It is known that G m contains m 3 
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|0)(0| 
|0)(0| 




Figure 1: Quantum circuit for the implementation of the dihedral POVM. 



elements [34]. POVMs that are covariant with respect to the Weyl-Heisenberg groups have 
a physical motivation. Since the position and momentum of a particle cannot be measured 
simultaneously by any projection- valued measurement one has to construct POVMs which 
measure both observables with a certain inaccuracy. This idea has already been described 
in [19]: starting from a wave packet, i.e., a unit vector 6 L 2 (M.) we define a set {M s j} 
of operators by 

M 8 , t := ^-e isP+ ^\^)(iP\e- isp - tQ . 

where s,t £ M and P and Q are the position and momentum operators, respectively. Ex- 
plicitly, they are defined by (Pip)(x) := —i(d/dx)ip(x) and (Qip)(x) := xip(x). We then 
have that 

/ M s jdsdt = l. 

J s,t 

The POVM {M Si t} provides an approximative realization of the classical phase space since 
the measurement outcome (s, t) can be interpreted as the point (s, t) in the phase space. In 
the following we are interested in finite dimensional approximations of this. Assume that we 
want to measure the position and crystal momentum of a particle on a lattice with m points 
for m = 2 k [35]. Furthermore, we assume that it is possible to transfer the state of such a 
system into k qubits of a quantum register. That means that we can implement a bijection 
of the basis states with Hamming weight one to the basis states of the Hilbert space C m 
of the k qubits. The canonical basis states \j) of C m denote the position eigenstates. The 
states corresponding to the state vectors Yl^o e 27Ttl ^ m \j) with I = 0, ... ,m — 1 are the 
eigenstates of the crystal momentum. Explicitly, the crystal momentum p can be defined by 
p := 2Trl/m — -iT. With this definition the values of p are in the interval [— tt, it] that meets the 
usual physical intuition of the one-dimensional Brillouin zone of an infinite one-dimensional 
crystal. Here we characterize the position and momentum simply by the integer values 
j, I = 0, . . . , m — 1. The cyclic translation of the position is given by the action of S m and 
a change of crystal momentum by the action of T m . Consider a rank-one positive operator 
with the property that neither the position nor the momentum of the corresponding 
state is completely undefined. Set 

M U :=^SLT l m ^){MT m l S m K 
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Due to irreducible group action the equation . l Mjj = l m holds and the operators 
Mjj define a POVM. For large m we can find states with corresponding state vectors 

such that both values j and I are approximately defined. Here the word "approxi- 
mately" is understood with respect to the cyclic topology, i.e., m — 1 and are "almost" 
the same value. A good choice for the POVM will be the following. Set := Ylj c j\j) 
where the coefficients Cj are chosen such that the function j i— ► \cj\ 2 has a unique max- 
imum at jo and the modulus of the values Cj decrease with increasing distance from jo 
in the cyclic topology. If all values Cj are real and they decrease not too quickly the 
momentum I of the state is around jo, too. Then the measurement values j, I can di- 
rectly be interpreted as a good estimation for the position and momentum values. We 
will show that an efficient implementation of the POVM can be found in the case where 

= l/^(l, Q!) a 2 ,... ) Q ! m / 2 - 2 ,a m / 2 - 1 ,Q ! m / 2 - 1 ,a m / 2 - 2 ,... ! a 2 ,a,l) T G C m with a G C 
and an appropriate normalization factor l/y/R. 

In the following we consider the group-generated POVMs with respect to G m and the 
natural representation ip defined by <p(g) = g for all g G G m . This representation is irre- 
ducible. Therefore, following AlgorithmElfrom Section^we can set U = l m since l m decom- 
poses (p into a direct sum of irreducible representations. The vector = (v\, . . . ,v m ) T G 
C m with the normalization |f i| 2 + . • • + |f m | 2 = 1/m leads to the POVM where the defining 
matrix M G C mxn is given by 

( V\ Vl ... Vl v m ... v m ... v 2 ... v 2 \ 

IT) — 1 77") — 1 >yri — 1 

t>2 V2U1 ■ ■ ■ Vl<jJ Vl ... Vi<jJ ... V3 ... V3U) 

\ V m V m U m ~ l . . . V m UJ V m -\ . . . V m -lUJ ... Vl ... ViU J 

Note that we identify vectors g\^) and h\^) for different g,h G G m that are equal up 
to a global phase factor. Consequently, the POVM consists of at most n = m? different 
operators. For example when m = 4 the vector \^>) = (vi, v%, V3, v^) T G C 4 with \vi\ 2 + 
|^2 1 2 + 1^3 1 2 + I ^4 1 2 = 1/4 leads to the POVM with n = 16 operators and the corresponding 
matrix M G C 4x16 where M is defined by 

/ V\ Vl Vl Vl Vi Vi V4 . . . V2 V2 t>2 V2 \ 

V2 V2% —V2 —V2% Vl V\i —Vl —V\% ... V3 V3I —V3 —V3I 
V3 -V3 V 3 -V3 V 2 —Vi V 2 -V2 ... V4 -V4 V4 -V4 

\ Vi —v^i —V4 V4.1 V3 —V31 —V3 V31 ... vi —v\i —vi v\i J 

The symmetry of M G C mx ™ can be described on the generators by the equations T m M = 
M(l m <S> S m ) and S m M = M(S m <S> T m ). Therefore the representation tp mon : G m — > U{n) 
is defined by (p m0 n(T m ) = l m (8> S m and (p m0 n(S m ) = S m ® T m . The symmetry of M can 
also be written as 

T m M = M(l m ® T m ) Fm8F '" and S m M = M(Tl S m ) F ™® F ™ 
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where we use the notation A x = XAX^ and the Fourier transform F m as defined in Section 
1(5. II We can write (l m <g> T m ) and (T m ® S m ) as direct sums 



T m ) — T m © T n 



T m and (T m ®S m ) = S m ®u m - 1 S ri 



ujS n 



By using the equations T m S m T m = coS m and (l m (g> S' m ) z = (Tm © S 1 ™) we can conjugate 



these matrices with the diagonal matrix Z = l m © T™ 1 



z _ c T t 



T„, © T m in order 



to obtain the equations 

T m M = M{l m ®T m )^ F ^ z and S m M = M(l m ® S m )^ F ^ z . 

These equations show that we have the decomposition Wip mon W^ = tp © ... © tp with 
the matrix W = Z\F m © F m ). The representation WVmonW''' contains m components 
tp. Following Algorithm from Section ^ we have to find a representation (p' that leads 
to the direct sum tp®tp' = tp®...®tp with m components tp. Consequently, we choose 
tp' = tp © • • • © <~p with m — 1 components <p. We now consider the extension of the matrix 
C = MW j = M(F m © F m )Z G C mxn to a unitary matrix C G W(n). The matrix C is an 
element of the intertwining space 



Int((/?, ip 



ip) = . . . , a„) © l m : a,- G 



C 



More precisely, we have C = ((yjmvi, . . . , \fmv m )F m ) © l m G C mxn . For example, with 
m = 4 we have the group G4 = (Si,T4) with S4 = (1,2,3,4) and T4 = diag(l, i, — 1, — i) 
that contains 64 elements. In this example we have the equation 



C 



(vi,v 2 ,v 3 ,v 4 ) 



\ 



( 1 

1 
1 

V 1 



1\ 

-1 



© 



/ 1 



V 



\ 



1 / 



G C 



4x16 



The matrix C G C mxn determines the first m rows of the matrix C we want to construct. 
The matrix C is a unitary matrix of the intertwining space 



lnt(p © . . . © tp, ip © . . . © ip) = {A © l m : A G 



l ] C C 



When we write C = ^4 © l m then the matrix C determines the first row of A. Explicitly, 
the first row of A is 

(y/mvi, \frnvm) F m . (3) 
The operation for the implementation of the POVM is defined by 

Mt = w*&(U © V) = (F m © i^)Z(At © l m ) G W(n). 

In this equation we lictve — l n _ m leading to U © V = l m © l n _ m — l n . The general 
scheme for the implementation of the matrix M' is shown in Figure [2j For m = 2 k the 
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|0)<0| 
|0)<0| 
|0>{0| 
|0>(0| 



- F m 



n / : rp—1 rp—2 rp—rn/A rp—m/2 T-ij" 

F * ' - L m m "' tyi Ttx "m 



{A ] ® I rn ) Z (F m ® F^) 

Figure 2: Circuit for the implementation of the POVM with respect to the Weyl-Heisenberg 
group and the vector \^) = (vi, . . . ,v m ) T . The vector |^) determines the matrix . 



circuit contains the k controlled operations 

rp—l rp—2 rp — m/4 rp—m/2 

for the implementation of the matrix Z. The matrix T m = diag(l, u>, w 2 , . . . , w m_1 ) can be 
written as Kronecker product 

r TO =(J JU )«-«(£ °)eW(m). 

Therefore, the matrices T m of the circuit in Figure [21 can be implemented efficiently on a 
register of qubits. 

The circuit in Figured is efficient if the matrix A that contains the vector © as first 
row can be implemented efficiently. We can find such a matrix for the POVM with the 
vector 

\m = _L (1, a a 2 . . . a ™/2-2 Q m/2-l a m/2-l ^m/2-2 ^2 Q -gT £ c? n (4) 

where we have a G C and the normalization k = 2m(l + |a| 2 + |a| 4 + . . . + |a| m_2 ). A 
matrix A £ U(m) that contains the vector © as first row is given by 

A = J L/2 i B ra/i ® 5 m/8 <g> . . . B 4 ® £ 2 ® Si (g> So) J m / 2 -^ 
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Figure 3: Implementation of the matrix A' where A is a matrix that contains the vector 
(J3J) as first row. This matrix is part of the circuit in Figured for the vectors (J3J. 



where we use the unitary matrices 



B: 



1 a> 
-1 



eW(2). 



Here is defined to be the permutation matrix which maps 2i \— ► % and (2i — 1) 
i = 0, . . . , k. In our example with m = 4we have the matrix 



-i for 



4 (B l ® B ) J 2 



x/2 + 2|a| 



/ 1 


a a 


1 


\ 


a 


-1 -1 


a 




a 


-1 1 


—a 




V i 


a —a 
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The circuit scheme for the implementation of the matrix 



At = F m A 



B 



m/2 \ m/4: 



B ] m/8 ® . . . ® B[® B ] 2 ® B[® Bq \ J m/2 



is shown in Figure |SJ 



7 Conclusions and outlook 

We have shown that a group-covariant POVM can be reduced to an orthogonal measure- 
ments by a unitary transform which is symmetric in the sense that it intertwines two dif- 
ferent group representations. The symmetry of the unitary transform can be used to derive 
decompositions which in several cases of interest (as the Heisenberg-Weyl group) leads to 
an efficient quantum circuit for the implementation of the POVM. 
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We have argued that POVMs are often necessary in order to understand why large 
quantum systems show typically classical behavior on the phenomenological level. The 
POVM with Heisenberg-Weyl symmetry as well as the example in [5] show that the POVMs 
which appear in this context are often covariant with respect to some group. 

Besides the physical motivation to study implementations of POVMs by means of or- 
thogonal measurements in terms of quantum circuits there is also a motivation from com- 
puter science. The so-called hidden subgroup problem [36] is an attractive generalization of 
the quantum algorithms for discrete logarithms and factoring [37]. The standard approach 
for the hidden subgroup problem consists in a Fourier transform for the respective group 
followed by a suitable post-processing on the Fourier coefficients [38]. For abelian groups 
this post-processing consists simply in an orthogonal measurement in the computational 
basis. However, for non-abelian group measurements which are in fact POVMs are often 
more advantageous, sec c. g. [39]. The POVMs which appear to be useful to solve hidden 
subgroup problems for non-abelian groups are naturally group-covariant. The methods pre- 
sented in this paper might be useful to find quantum algorithms for the hidden subgroup 
problem for new classes of non-abelian groups. 
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